Demonstration of reversible phase-insensitive optical amplifier 
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We experimentally demonstrate phase-insensitive linear optical amplification which preserves the 
idler at the output. Since our amplification operation is unitary up to small excess noise, it is 
reversible beyond the classical limit. The entanglement between the two output modes is the resource 
for the reversibility. The amplification gain of G = 2.0 is demonstrated. In addition, combining this 
amplifier with a beamsplitter, we also demonstrate approximate cloning of coherent states where 
an anticlone is present. We investigate the reversibility by reconstructing the initial state from 
the output correlations, and the results are slightly beyond the cloning limit. Furthermore, full 
characterization of the amplifier and doner is given by using coherent states with several different 
mean values as inputs. Our amplifier is based on linear optics, offline preparation of nonclassical 
ancillas, and homodyne measurements followed by feedforward. Squeezed states are used as the 
ancillas, and nonlinear optical effects are exploited only for their generation. The ancillas introduce 
nonclassicality into the amplifying operation, making entanglement at the output. 

PACS numbers: 03.67.Hk, 42.50.Dv, 42.65.Yj 
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I. INTRODUCTION 

Quantum optics is governed by rules imposed by com- 
mutation relations which have to be kept during time 
evolution. Optical amplification is no exception to this 
story. Typically, the amplified output is suffered from 
inevitable excess noise. This limitation is quantum- 
mechanically imposed, thus does not depend on the spe- 
cific realization methods. Caves classified general linear 
amplification into phase-insensitive amplification (PIA) 
and phase-sensitive amplification (PSA) [l| . He also sys- 
tematically derived the quantum limit of excess noise for 
such general linear amplification with arbitrary gain from 
the requirement to preserve commutation relations. This 
excess noise originates from quantum fluctuations in the 
auxiliary system required to keep energy conservation. 

We concentrate on PIA, supposing the target of am- 
plification to be optical wave amplitude of a single mode, 
which is denoted by the term "signal". Classical coun- 
terpart of PIA is a conversion of arbitrary complex wave 
amplitude a e C into \[Gu, where G > 1 is the gain of 
amplification. As is found in ordinary textbooks, anni- 
hilation operators in quantum optics correspond to com- 
plex amplitudes in classical optics. Therefore, we de- 
scribe the amplifying process by the transformation of 
annihilation operators. Quantum-mechanically optimal 
PIA in the sense that the excess noise is minimized can 
be achieved by the following transformation [l[: 



(i) 



where a"? and 5°^* are the signal mode's annihilation op- 
erators before and after the amplification, respectively. 



There is an extra term e %ti \JG — 1 (a^i)^ which is in- 
troduced in order to meet the commutation relation of 
[a s j g , a sig ] 
modes. Here, 



1 for both the input and output signal 
G 1 is an arbitrary phase factor, and 



a|^[ is another mode's annihilation operator in the aux- 
iliary system. Throughout this paper, the ancilla mode 
represented by a-gj is denoted by the term "idler" and 
distinguished from other ancilla modes. Eq. ([T]) becomes 
the input-output relation of optimal PIA when the idler 
is in a vacuum state. The quantum fluctuation of the 
idler contaminates the amplified signal. This is the in- 
evitable excess noise of PIA. Note that this penalty pre- 
vents amplification from being a loophole of the uncer- 
tainty relation in joint measurements @, Q. At the limit 
of high-gain amplification, we can see the famous 3 dB 
cost of the noise figure for PIA of coherent states. In 
addition to this intrinsic excess noise, further nonintrin- 
sic excess noise may be caused by other ancilla modes in 
nonoptimal PIA. 

There are numerous practical realizations of optical 
amplification. Doped fiber amplifiers (DFAs) and semi- 
conductor optical amplifiers (SOAs) utilize stimulated 
emissions [J], and Raman amplifiers (RAs) and optical 
parametric amplifiers (OPAs) utilize nonlinear optical ef- 
fects. In principle, there is no quantum-mechanical rea- 
son to prevent these realizations from achieving the opti- 
mal PIA in the form of Eq. ([1]). However, the real devices 
with current technology are accompanied by further ex- 
cess noises. 

Recently, PIA operating almost at the optimal level is 
experimentally demonstrated by Josse et al. by utiliz- 
ing feedforward 5] . The reason for the high efficiency of 
Josse's PIA is that it does not require inefficient nonclas- 
sical operations or nonclassical ancillas. It uses a vacuum 
state as an ancilla which is present everywhere, and linear 
optics and homodyne measurements followed by feedfor- 
ward which are highly efficient. 

Although Josse's PIA is a good attainment, it is not 
the end of the story. The signal transformation in Eq. (fTJ) 
is an irreversible thermalizing process. Complete PIA 
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should have unitary realization on an expanded Hilbert 
space. In order to unitarize PIA, two-mode description 
is sufficient. The full input-output relation becomes as 
follows: 

a°£ =VGa% + e ie VG=l&, (2a) 
agf =VGa* + e* e VG^T(a%)l (2b) 

Note that the roles of the signal and idler are symmetric 
in this relation. 

The significance of unitarization must be the reversibil- 
ity. The inverse transformation is easily derived when we 
take notice of the fact that Eq. ([2]) is equivalent to two- 
mode squeezing operation. A two-mode squeezing oper- 
ation parametrized by (G, 9) is canceled by another two- 
mode squeezing operation where the squeezing direction 
is opposite, i.e., (G, 9 + ir). Nonetheless, in many ampli- 
fication schemes including Josse's experimental demon- 
stration Q, the idler output is lost in the inextractable 
environment, making the process irreversible. 

In order to realize idler-preserving and close-to-optimal 
PIA, we require some nonclassicality for the amplifier. 
This is contrastive to Josse's idler-nonpreserving PIA 
which does not require any nonclassicality. A typical 
strategy to introduce nonclassicality into feedforward- 
based quantum circuits is to use nonclassical states as 
ancillas. Continuous-variable (CV) quantum teleporta- 
tion [6| and CV error correction [7( are good examples. 
In these examples, squeezed states are used as ancillas 
that support the performance beyond the classical limit, 
and the complex operations after the state preparation 
stage arc efficiently implemented by linear optics. 

In this paper, by employing the feedforward-based 
scheme proposed in Ref. @], we demonstrate PIA of 
coherent states which preserve the idler output. The 
scheme basically relies on linear optics including homo- 
dyne measurements and feedforward. Squeezed vacuum 
states are used as ancillas, which inject nonclassicality 
into our PIA. Only for generating the nonclassical an- 
cilla states, we resort to nonlinear optical effects. Our 
demonstration is for the amplification gain of G = 2.0, 
which is tuned via passive optical devices and feedforward 
electric circuits. Combining PIA for G — 2.0 with a half 
beamsplitter, we also demonstrate 1 — > 2 approximate 
cloning of coherent states, where an "anticlone" remains 
at the output. (Anticlone will be explained in Sec. Hill ) 
In principle, our amplifier and doner becomes quantum- 
mechanically optimum at the limit of infinite squeezing 
of the ancillas. For the case of finite squeezing, as is the 
real situation in experiments, further excess noise invades 
in accordance to the level of the squeezing. However, the 
degradation is small enough to retain nonclassical fea- 
tures. The behaviors of our amplifier and doner are fully 
characterized by using several coherent states as inputs. 
Furthermore, we also pay much attention to the output 
correlations, because nonclassical properties clearly ap- 
pear in them. For the PIA experiment, we check the 
Einstein-Podolsky-Rosen (EPR) correlation between the 



signal and idler outputs. For the cloning experiment, we 
check bipartite entanglement between each clone and the 
anticlone, which as a whole proves tripartite entangle- 
ment of class I 0. Moreover, for both experiments, the 
reversibility is investigated from the output correlations. 

Our idler-preserving PIA is significant in several re- 
spects. First of all, the reversibility will pave the way 
to new schemes. Recently, there is a proposal of a CV 
quantum interface that enables in principle a unit fidelity 
of transfer using such reversible PIA Moreover, the 
reversibility in cloning is also advantageous. Cloning of 
unknown states is distribution of information, and its re- 
versibility reserves the option to recover the distributed 
fragments of the information. This will be further dis- 
cussed in Sec. Hill Secondly, our PIA would have some 
applications as two-mode squeezing operation. Note that 
one-mode squeezing operation is already demonstrated 
successfully in Ref. [ll| with similar approach. 

In this introduction, PIA has been described together 
with a brief historical review. Especially, the nonclassical 
property of PIA is discussed, which is obscure in many 
amplification processes because the idler output is lost in 
the inextractable environment. The subsequent contents 
of this paper are as follows. In Sec. [HI feedforward-based 
PIA is described, explicitly showing the excess noise due 
to finite squeezing of ancillas. In Sec. Hill CV quantum 
state cloning and its connection with PIA are described. 
In Sec. HVl the experimental setup is described. In SeclVl 
the experimental results for PIA of coherent states with 
G = 2.0 are shown. In Sec. IVI1 the experimental results 
for 1 — > 2 approximate cloning of coherent states are 
shown. In Sec. IVII1 our experimental achievements are 
summarized. 



II. FEEDFORWARD-BASED AMPLIFIER 

In our definition, feedforward means that the opera- 
tions after some measurements are changed depending 
on the measurement outcomes which in general are ob- 
tained randomly. In particular, in this paper it indicates 
phase space displacement operations whose amounts are 
proportional to the results of homodyne measurements. 

We know two specific schemes for feedforward-based 
PIA that preserves the idler at the output. One scheme 
is proposed by Filip et al. in Ref. Q, in which PIA is 
composed of two feedforward-based single-mode squeez- 
ers proposed in the same paper. The other scheme is 
proposed by Josse et al. in Ref. [| as a modification 
of the idler-nonpreserving PIA. Note that Josse's idler- 
preserving PIA is just a theoretical proposal and the 
idler-nonpreserving PIA alone is experimentally demon- 
strated. 

Both of Filip's scheme and Josse's scheme rely on lin- 
ear optics including homodyne measurements and feed- 
forward, and require offline-prepared nonclassical states 
as ancillas. Moreover, in both schemes, the gain of 
amplification is accurately and stably determined via 
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the choice of passive optical devices and correspondingly 
feedforward gains. As for the nonclassical ancillas, Filip's 
scheme requires two single-mode squeezed states, on the 
other hand, Josse's scheme requires a two-mode squeezed 
state. Since two single-mode squeezed states can be con- 
verted to a two-mode squeezed state and vice versa by 
a half beamsplitter interaction, the amounts of the non- 
classical resources required for the two distinct schemes 
are the same. 

For both schemes, the feedforward-based PIA coin- 
cides with the quantum-mechanically optimal PIA only 
at the limit of infinite squeezing of the ancillas. For the 
case of finite squeezing, excess noise contaminates the 
output to some extent. Note that this is a common 
matter of feedforward-based CV deterministic process- 
ing 0, @, [Hi . The difference between the two schemes 
proposed by Filip and Josse solely arises in this excess 
noise. For Filip's scheme, it appears symmetrically in the 
signal and idler outputs. On the other hand, for Josse's 
scheme, it appears only in the idler output. The better 
choice between different schemes depends on the specific 
application. 

We have chosen the symmetrical one. In the demon- 
stration in Sec. [Vj we confirm the symmetry of PIA by 
swapping the roles of the signal and idler. 



Fig. 1(a) shows the schematic of our PIA, from which 
the symmetry of the signal and idler is obvious. Its de- 
tails will be described in Sec. lIVI Here we give the input- 
output relation. In the following, the quadrature phase 
amplitudes of each optical mode are denoted by x and 
p, which correspond to the real and imaginary parts of 
the mode's annihilation operator a, i.e., a = x + ip. The 
phase factor 8 in Eq. ^ can be arbitrarily changed by 
pre- and post-processing of phase rotation of the idler. 
Therefore, we consider the case of 8 = without loss of 
generality. Explicitly showing the excess noise coming 
from finitely squeezed ancillas, the input-output relation 
becomes as follows |8|: 



=M^ + ^)^i n + M^-^)4 n -V^A t , (3a) 
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The subscripts '1' and '2' represent the two main modes. 
They correspond to the signal and idler, though we do 
not specify which is which because the relation is sym- 
metric. £a and pb denote the squeezed quadratures of 
the two ancilla modes. At the limit of infinite squeezing, 
these terms vanish, and the transformation above strictly 
coincides with the optimal PIA. The amplification gain 
G is determined via one parameter i?, which is a com- 
mon reflectivity of two beamsplitters BS-A and BS-B in 
Fig. l(a)| with the relation of 

G =i(7I + ^) 2 - ( 4 ) 



One-to-one correspondence of 1 < G < oo and < R < 1 
is easily checked. Note that the feedforward gain is also 
parameterized by R. It is chosen so that the antisqueezed 
noises from the ancillas are canceled out at the output. 

For the demonstration of G = 2, The value of R should 
be chosen as 3 — 2^/2 « 0.17. The resulting input-output 
relation becomes as follows: 



=v / 2 4 n + x I 1 n + 



-out 
X 1 



V2- iil ut , 



p° 2 ut =V2p^- PT 



V2-lpg\ 
V2-lxT\ 
V2~lp° B ut . 



(5a) 
(5b) 
(5c) 
(5d) 



III. QUANTUM STATE CLONING 

It is known as the no-cloning theorem that an un- 
known quantum state 1^) cannot be perfectly duplicated 
as l^)!^) [13j |. However, approximate cloning is possible, 
which can go beyond some classical limit in general. 

In this section, we will discuss CV cloning, and make 
its connection with PIA. Furthermore, its reversibility 
is discussed by introducing the notion of anticlone. In 
general, cloning can be described as a unitary operation 
supported by ancilla systems. The ancilla output sys- 
tem generally depends on the cloned state, from which 
anticlones are obtained. We show the equations for 
1—7-2 cloning, which corresponds to the experimental 
demonstration in Sec. IVII However, PIA allows general 
K — > L cloning in principle, which will be described in 
Appendix [A] Here, the notation K — > L means that L 
clones are created from K identical originals. 

First, we would like to say that pragmatic cloning for 
CV is not the universal cloning |14j with respect to the 
infinite-dimensional Hilbert space, because it is an unnat- 
ural situation that all the states in the noncompact space 
appears with equal probability. In general, the choice of 
the appropriate doner depends on how the information 
is embedded in the infinite-dimensional Hilbert space. 

The typical situation is that the CV information is 
embedded as a displacement on some quantum state 
Here, which we refer to as a core state, is either 
known or unknown. Then, the set of possible origi- 
nal states is S = {D(xd,Pd)\4>) I ( x d,Pd) £ K 2 }, where 
D(xd,Pd) = exp[— 2i(xdp — Pax)] is the displacement op- 
erator. [The probability density p(xd,Pd) is omitted be- 
cause we consider for simplicity the case where (xd,Pd) is 
uniformly distributed.] As a special case of this, the set 
S becomes all coherent states when the core state \ip) is 
known to be a vacuum state. This way of embedding is 
found in ordinary CV quantum key distribution (QKD) 
protocols [lj|. 

For such protocols, the role of cloning is distribution 
of the information, rather than duplication of a quantum 
state. Therefore, the measure of the cloning precision 
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should be related to the estimation of {xd,Pd)i instead of 
the traditional fidelity. Furthermore, asymmetric doner 
is significant as well as symmetric doner. Arbitrary share 
ratio of the information is achieved by doner with tunable 
asymmetry. 

We suppose a simple picture of cloning where some 
noise is added to the original state as the penalty of 
cloning. Then, the quality of cloning is totally deter- 
mined by this noise. For simplicity, we impose rotational 
symmetry on the noise added to each clone. This is nat- 
urally justified when the core state \tp) is either known 
to be symmetric or unknown. The added noise is charac- 
terized by its variance n k = (Aa££f£) 2 + ( A Pcin-f ) 2 0» 
where k £ {1, 2} for 1 — > 2 cloning. Note that n k corre- 
sponds to the mean photon number of thermalization in 
the fc-th clone. 

The variances Ufc are directly connected to the mean 
square errors in the estimation of (xd,Pd)- Therefore, a 
measure can be constructed from them. Given the de- 
sired asymmetry, the cost function is determined (l6j : 

C(ni,n 2 ) = c\n\ + c 2 n 2 . (6) 

The positive parameters Ck determine the asymmetry, 
(ci = C2 corresponds to symmetric cloning.) The doner 
that minimize the cost function is the optimal. 

It is obvious that the optimal doner is Gaussian when 
the cost function is set as a function of the noise vari- 
ances as in Eq. ([5]). In order to minimize it, the ancillas 
that support cloning are chosen in minimum uncertainty 
states, which are Gaussian. This contrasts with the eval- 
uation by the fidelity. Non-Gaussian cloning can slig htly 
go beyond the Gaussian fidelity for coherent states [171 ]. 
We emphasize that the evaluation by the variances is 
more practical. It is pointed out that optimal attack in 
QKD is Gaussian [ljj. 

There is a restriction on the excess noises ilk which is 
imposed by quantum mechanics pol l2lj : 

mnj > (1/2) 2 . (7) 

The optimal doner with respect to the cost function in 
Eq. ([5]) satisfies the equality in Eq. (J7J by necessity. This 
noise penalty comes from consistency with the uncer- 
tainty relation. The attainable information of the origi- 
nal state does not increase by cloning due to this noise. 
Recall that the inevitable noise in PIA comes from the 
same reason. 

Indeed, the optimal doner can be constructed from 
the optimal phase-insensitive amplifier and beamsplit- 
ters. For example, 1 — > 2 cloning with arbitrary asym- 
metry is achieved by putting an amplifier in one of the 
arms of a Mach-Zehnder interferometer [2l|. Especially, 
for symmetric cloning, the reflectivity of the first beam- 
splitter becomes unity, i.e., it is achieved by first am- 
plifying the original with G = 2 and then splitting the 
amplified signal in half. This procedure can be extended 
to K —> L cloning [l6|, |22| . The optimality of this re- 
alization is proven with respect to the cost function in 
Eq. © 0. 



For Gaussian cloning of coherent states, the added 
noise variance rik and the fidelity have correspondence 
as Fk = 1/(1 + rik). By using Eq. 0, the upper limit of 
fidelity is obtained for arbitrarily asymmetric Gaussian 
cloning. In particular, it becomes F = 2/3 for the sym- 
metric case. This is significantly higher than the classical 
limit of F — 1/2, where we regard the limit of state esti- 
mation as the classical limit of symmetric cloning because 
the estimated state is classical information which can be 
copied any number of times. Note that the sameness of 
state estimation and asymptotic cloning where the num- 
ber of clones tends to infinity is proven for a general set S 
of possible original states |23| . We refer to these fidelities 
only for the consistency with previous works. We stress 
again that our actual interest is the variances. 

We have seen above that the clones are made of the 
signal output of the amplifier. When cloning is unitarily 
realized, we still have the idler output, whose state is 
affected by the original state. Now we pay attention to 
this ancilla output system. 

As mentioned at the beginning of this section, anti- 
clones are byproducts of cloning which are obtained from 
the ancilla systems. Especially, for 1 — > 2 cloning, the 
idler output itself is an anticlone. It is an approximation 
of the phase-conjugated original state, or in other words, 
the output of an approximate NOT gate. Qubit version 
of this gate is demonstrated in Ref. [2J| . 

Anticlones are important when we are concerned with 
the reversibility. The originals can in principle be per- 
fectly reproduced only when all the clones and anticlones 
are present. For the reversibility, the essential resource 
is nonclassical correlations, or entanglement. Concep- 
tually, the excess noises in cloning is canceled by using 
the nonclassical correlations. Therefore, we discuss the 
existence of entanglement in the three-mode output sys- 
tem of 1 — > 2 cloning. Clones are obtained by split- 
ting the amplified signal, thus there is no entanglement 
among clones. However, there is entanglement between 
each clone and the anticlone. We stress that the resource 
for the recovery is not the anticlones themselves but the 
entanglement. We can obtain anticlones without entan- 
glement with clones as follows. Suppose the situation 
where two independent doners are running, and the same 
states are used as the inputs for them. Then, the clones 
obtained from one doner do not have entanglement with 
the anticlones from the other. In this case, the originals 
can not be recovered from the noncorrelated outputs. 

For the recovery of the originals, the inverse unitary 
operation is not required. The optimal cloning can be 
fully reversed by the Bell measurement on a clone and 
an anticlone and subsequent feedforward to the remain- 
ing single clone (25[. Note that this recovery scheme 
works not only for coherent states but also for an ar- 
bitrary core state \ip). This scheme is efficient from two 
points of view. One is on a technical level that the homo- 
dyne measurements and feedforward displacement oper- 
ations are quite efficient with current technology. The 
other is on a conceptual level that the performer of the 
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Bell measurement and the owner of the remaining clone 
who is willing to recover the original can be spatially 
separated. For this case, they only need classical chan- 
nels for communication, and never quantum channels. 
Note that even partial reversal is possible with a similar 
scheme based on local operations and classical commu- 
nication (LOCC), which converts, e.g., symmetric clones 
to asymmetric clones [25j . 

We would like to discuss practical aspects of cloning 
and its reversibility assisted by classical communication. 
As described above, cloning of a quantum state is re- 
garded as distribution of information among plural par- 
ticipants. The information of the original is to some ex- 
tend accessible to individual participants. This situation 
is clearly distinguished from that found in usual quantum 
error correcting codes where the quantum information is 
mapped on a larger Hilbert space so that no information 
about the original is accessible from a localized system. 
Such share of information would play important roles in 
several scenarios, in which the reversibility would give a 
tactical aspect to information exchange. For example, 
cloning is a possible attack by an eavesdropper in QKD. 
In this example, the reversibility of cloning provides the 
opportunity for the communicators to negotiate with the 
eavesdropper when they know the attack (25[. Since co- 
herent states are a strong candidate for the information 
carrier in quantum communication, cloning of coherent 
states is especially of great significance. 

There are several experimental previous works which 
demonstrate cloning of coherent states beyond the clas- 
sical limit of F = 1/2 in non- reversible ways, i.e., their 
anticlones are lost in the environment. In Ref. [261 ] . us- 
ing feedforward-based PIA of Ref. [5| , almost quantum- 
limited 1 — > 2 cloning is demonstrated. In Ref. [27j . 
telecloning is demonstrated where the original coherent 
state is teleported and cloned at the same time. 

In Sec. lVIl we demonstrate 1 — > 2 symmetric Gaussian 
doner which preserves an anticlone at the output. As 
we apply a half beamsplitting to 
ihe feedforward-based PIA with the 
. [TTJ In the demonstration, 



is shown in Fig. 1(b) 



the signal output of t 
gain G = 2 described in Sec. 
the reversibility is checked from the output correlations. 
To our knowledge, there is no previous experiment of this 
kind even in qubit regime. Although our demonstration 
is only for coherent states, our doner should equally work 
for arbitrary core state \ip) as discussed above. 

We close this section by giving the input-output rela- 
tion of the optimal 1 — > 2 symmetric cloning. By substi- 
tuting G = 2 and 9 = into the input-output relation in 
Eq. ([2]) and splitting the signal output in half, we obtain, 



^cln-1 ^org T" idl ^~ "^/2 ac ' 
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(8a) 
(8b) 
(8c) 



respectively. The subscript 'idl' denotes the idler input 
for PIA which is in a vacuum state. The annihilation op- 
erators with the subscript 'vac' indicate another ancilla 
in a vacuum state, which invades from the empty port 
of the final half beamsplitter. For the excess noise of the 
two clones, n\ = = 1/2 is easily checked. Therefore, 
this doner is optimum when evaluated by the cost func- 
tion in Eq. ([5]) with c\ = C2- When PIA is realized with a 
feedforward-based scheme as is found in Sec. EH further 
excess noise contaminates the output in accordance with 
the squeezing levels of the ancillas. 



IV. EXPERIMENTAL SETUP 



Schematic of the experimental setup for PIA is illus- 
trated in Fig. l(a)[ and that for approximate cloning is 
illustrated in Fig. 1(b) The light source is a Ti:sapphire 



laser, which has a continuous-wave single-mode output 
of 860 nm in wavelength and about 1.5 W in power. We 
treat the quantum states of narrow sidebands located at 
1.34 MHz apart from the optical carrier frequency. 

Two main beams that go from in-1 and in-2 to out-1 
and out-2 in Fig. 1(a) carry the quantum states which 
are targets of PI A. The setup has a form of a Mach- 
Zehnder interferometer that holds a single-mode squeezer 
(squeezer-A or squeezer-B) in each arm. This decompo- 
sition of unitary PIA into squeezers and beamsplitters is 
derived from the bosonic version of Bloch-Messiah reduc- 
tion shown in Ref. [28| . We note that this setup is almost 
the same as that for quantum nondemolition (QND) in- 
teraction demonstrated in Ref. [2!|. This fact shows the 
capability of our setup to realize many types of two-mode 
Gaussian interaction. Combining PIA for G = 2.0 with 
another half beamsplitter as is shown in Fig. 1(b) 1—^2 



where the subscripts 'org', 'cln-1', 'cln-2', and 'a-cln' de- 
note the original, first clone, second clone, and anticlone, 



approximate cloning of coherent states is achieved. 

Squeezer-A and squeezer-B are feedforward-based 
squeezers, which are theoretically proposed in Ref. Q and 
experimentally demonstrated in Ref. [11]. Each squeezer 
consumes an ancilla in a squeezed state, which is gener- 
ated by an optical parametric oscillator (OPO). 

Note that several essential optical elements are omit- 
ted from Fig. [TJ such as a second harmonic generation 
(SHG) cavity to generate pump beams for OPOs, and 
three spatial- mode cleaning cavities (MCCs). One MCC 
is used for local oscillators (LOs) for homodyne mea- 
surements and auxiliary beams for feedforward displace- 
ments. The other two MCCs are used for individual input 
beams. 

The experimental procedure is divided into three steps: 
Firstly, we prepare input coherent states and ancilla 
squeezed vacuum states. Secondly, we implement PIA 
and cloning via feedforward. Finally, the output states 
are homodyne measured for verification. In the following, 
we describe the experimental details of each step. 
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(a)Sctup for PIA. 



(b)Setup for 1—^2 doner. 



FIG. 1: Schematics of experimental setups. 0P0: Optical parametric oscillator. EOM: Electro-optic modulator. LO: Local 
oscillator. HBS-P, HBS-F, HBS-C: Half beamsplitter. BS-A, BS-B: Beamplistter with reflectivity R « 0.17. Four beamsplitters 
(HBS-P, HBS-F, BS-A, and BS-B) are variable, composed of two polarization beamsplitters and one half-wave plate. 



A. Preparation 

At this step, we generate coherent states which arc 
used as inputs, and squeezed vacuum states which are 
used as ancillas. 

The nonzero mean values of the sideband coherent 
states at 1.34 MHz are produced by appropriately mod- 
ulating the optical carriers. In our setup, the relative 
phase of interference at each beamsplitter is designed to 
be fixed with active feedback control. Therefore, in order 
to make an arbitrary phase space displacement in the in- 
put modes, both amplitude modulation (AM) and phase 
modulation (PM) are utilized. AM and PM make non- 
zero mean values of and for the first input mode, 
and those of p™ and x™ for the second input mode, re- 
spectively. Each of four electro-optic modulators (EOMs) 
before PIA in Fig. 



1(a) corresponds to one of these four 



quadratures. On the other hand, in Fig. 1(b) only two 
EOMs are depicted before PIA which are both located at 
the first input beam path. Therefore, the symmetry of 
the two input modes is broken in the cloning experiment. 
One input mode is the target of cloning, and the other 
input mode is set in a vacuum state throughout. For 
both experiments, the modulations are switched on and 
off in order to use several coherent states as inputs. After 
these EOMs, there are MCCs though they are omitted 
from Fig. [T] 

Squeezed vacuum states are each generated by an 
OPO which is driven below the threshold. Our OPO 
has a bow-tie shaped configuration with a round-trip 
length of about 500 mm. It contains a periodically- 
poled KTiOP0 4 (PPKTP) crystal as a nonlinear optical 
medium, which is commercially available from Raicol and 
has 10 mm length and 1 mm by 1 mm cross section. The 
experimental datails of our OPO squeezing are found in 
Ref. [30( | . The squeezing level with the pump of about 
100 mW is about —5 dB relative to the shot noise level 
at 1.34 MHz. The pump beams for the OPOs are the 



second harmonic of a fundamental beam, generated by 
a SHG cavity. Most of the Tksapphire laser output is 
sent to the SHG cavity, whose output of about 300 mW 
is divided into two to pump the individual OPOs. The 
SHG cavity has almost the same configuration as that of 
OPOs, whereas a KNb03 (KN) crystal is used instead of 
the PPKTP crystal. 

Modulation sidebands other than 1.34 MHz are ex- 
ploited for active feedback control of optical interferences. 
A modulation at 13.5 MHz is utilized for locking cavities, 
including the SHG cavity, the two OPOs, and the three 
MCCs. On the other hand, lower frequency modulations 
at 193 kHz and 333 kHz are utilized at the OPOs to lock 
the phases of the pump beams. Furthermore, the two in- 
put beams are modulated at 108 kHz and 154 kHz. These 
four low-frequency modulations contribute to the lock of 
the downstream interferometric system in the subsequent 
steps, as mentioned later. 



B. Amplifier and Cloner 

The two input beams are combined at a preceding 
half beamsplitter (HBS-P) and then sent to squeezer-A 
and squeezer-B. After the squeezing operations, the 
two beams interfere again at another half beamsplitter 
(HBS-F), which completes PIA. By splitting one of the 
two output beams by another half beamsplitter (HBS-C), 
l->2 cloner is obtained. 

The squeezing procedure goes as follows. First, the 
main beam is combined with an ancilla beam coming 
from an OPO at a beamsplitter (BS-A or BS-B). Next, 
one of the two beams after the beamsplitter is homo- 
dyne measured. Finally, the measurement outcome is 
fed forward to the remaining beam. BS-A and BS-B have 
the common reflectivity of R. This parameter R deter- 
mines the degree of the feedforward-based squeezing and 
thus the gain of amplification G with the relation shown 
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in Eq. (|4|. As is already mentioned, R ~ 0.17 for our 
demonstration of G = 2.0. 

The feedforward operation is a phase space displace- 
ment whose amount is proportional to the random out- 
come of the homodyne measurement. The electric sig- 
nal from the homodyne detector is sent to an EOM to 
be converted into an optical signal, where the gain and 
phase at 1.34 MHz are carefully chosen. The auxiliary 
beam which is modulated by the feedforward EOM has 
the power of 150 /uW, 1% of which subsequently enters 
the mainstream via an asymmetric beamsplitter (99:1). 

The powers of the two input beams are 10 fiW, and 
those of the two ancilla beams are 2 /iW. These pow- 
ers are considerably smaller than 3 mW of LOs used for 
homodyne detections. 

The four beamsplitters of PIA (HBS-P, HBS-F, BS-A, 
and BS-B) are actually composed of two polarization 
beam splitters and a half- wave plate in the same manner 
as the QND experiment in Ref. [29J. Their reflectivities 
are variable by rotating half-wave plates. They enable us 
to measure the input states as well as the output states 
with the same homodyne detectors for verification. The 
propagation losses of two main beams are measured to be 
7% on average, which mostly come from these variable 
beam splitters. 

In order to control the relative phases at beamsplitters 
with active feedback, interferences between the carriers 
and the low-frequency modulations are monitored. This 
is typically done by picking up 1% of the beam after the 
interference, though such details are omitted from Fig.[T] 
For each locking point, an appropriate modulation side- 
band is chosen, and the error signal is extracted from 
the interference between the carrier and the sideband 
by demodulation. However, HBS-P and HBS-F are ex- 
ceptions, where the interference between two modulation 
sidebands, namely 108 kHz and 154 kHz, is exploited. 
The beat frequency of 46 kHz is chosen for the reference 
signal of demodulation to obtain the error signals. 



C. Verification 

PIA is characterized by measuring two-mode input 
states as well as two-mode output states using two ho- 
modyne detections. In the cloning experiment, on the 
other hand, three-mode output states are compared with 
single-mode input states. The input states are measured 
by setting the reflectivities of the four variable beamsplit- 
ters to unity and disabling the feedforward. The quan- 
tum efficiency of a homodyne detector is about 99%, and 
the dark noise is about 17 dB below the optical shot noise 
produced by the LO. The interference visibilities to the 
LOs are 98% on average. 

The outcomes of the final homodyne measurements are 
analyzed in either of the two ways below. 

In one way of analysis, the quadrature data are di- 
rectly treated, which are obtained by lock-in detection 
of 1.34 MHz components of the homodyne outputs. A 



signal from a homodyne detector is mixed with the refer- 
ence signal at 1.34 MHz, and then low-pass filtered with 
the cutoff of 30 kHz. Subsequently it is analog-to-digital 
(A/D) converted for storage with the sampling rate of 
300 kHz and the resolution of 14 bits (PXI-5122, National 
Instruments Corporation). In this analysis, the phase of 
the homodyne detection is slowly scanned. The phase in- 
formation is stored simultaneously with the quadrature 
values using the same A/D board. From the resulting 
marginal distributions, phase space distributions (i.e., 
Wigner functions) are reconstructed, where we assume 
that all the quantum states obtained in the experiments 
are Gaussian. The first and second moments are com- 
puted so that the likelihoods are maximized. 

The other way is the power analysis at 1.34 MHz us- 
ing a spectrum analyzer. In this analysis, the measured 
quadratures are set to either x or p. Not only the powers 
of the output quadratures but also those of their correla- 
tions are measured for several input coherent states. The 
resolution bandwidth is 30 kHz, the video bandwidth is 
300 Hz, the sweep time is 0.1 s, and 20 times averaging 
is taken for each trace. 

Note that we can easily see the effect of the Hermitian 
conjugate term in Eq. (|2b[) as a mirror image with the 
former way of analysis, whereas we can not do this with 
the latter. 



V. EXPERIMENTAL RESULTS FOR 
PHASE-INSENSITIVE AMPLIFIER 

The two main modes are denoted by "mode-1" and 
"mode-2" . One of them is the "signal" and the other is 
the "idler" , which are initially in a coherent state and a 
vacuum state, respectively. By swapping the role of the 
signal and idler, we check the symmetry of our PIA. 

We first show the results of the lock-in detection, be- 
cause it is intuitively easier to see. Figs. [5] and 0] show 
the experimental quadrature values at various phases of 
LOs. For Fig.O the mode-1 is the signal and the mode-2 
is the idler, whereas for Fig. [2J the mode-2 is the sig- 
nal and the mode-1 is the idler. There are three sub- 
figures corresponging to the signal input (a), the signal 
output (b), and the idler output (c). Horizontal axes 
arc the measurement phases <f> which are scanned from 
to 2tt. The quadrature at <j) — corresponds to x and 
that at <j> = ir/2 does to p. Vertical axes are normal- 
ized quadrature values where the standard deviation of 
vacuum fluctuation is 0.5. Each set of data is taken for 
about 0.2 seconds. Quadrature data are plotted every 10 
points in the figures, whereas the whole data are used for 
the analysis. The sinusoidal curve of the signal input rep- 
resents the nonzero mean amplitude of a coherent state, 
and the fluctuation around the sinusoid represents the 
quantum noise. With regard to the fluctuation, it grows 
uniformly at both the signal and idler outputs. This uni- 
formity is an evidence of the phase-insensitivity of our 
amplifier. On the other hand, with regard to the sinu- 
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FIG. 2: Quadrature data. Mode-1 for signal and mode-2 for 
idler. The phases of homodyne measurements are scanned 
from to 27t, which correspond to horizontal axes. Vertical 
axes are quadrature values. 



FIG. 4: Quadrature data. Mode-2 for signal and mode-1 for 
idler. The phases of homodyne measurements are scanned 
from to 2-7T, which correspond to horizontal axes. Vertical 
axes are quadrature values. 
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FIG. 3: Phase space distributions computed from the quadra- 
ture data in Fig. [2] The first and second moments of Gaussian 
Wigner functions are represented by ellipses. Green: Signal 
input. Red: Signal output. Blue: Idler output. 



FIG. 5: Phase space distributions computed from the quadra- 
ture data in Fig. [4] The first and second moments of Gaussian 
Wigner functions are represented by ellipses. Green: Signal 
input. Red: Signal output. Blue: Idler output. 



soidal curve, the two output modes show different behav- 
iors. At the signal output, the amplitude of the sinusoid 
is amplified from that of the signal input, maintaining 
the phase. At the idler output, the amplitude of the si- 
nusoid is the same as the signal input, whereas the phase 
is flipped. This flip is due to the Hermitian conjugate 
term in Eq. (|2bl) . For both figures, the same qualitative 
behaviors as mentioned above are observed. 

Figs. [3] and [5] are phase space diagrams, which are com- 
puted from the quadrature data shown in Figs. [5] and 21 
respectively. The experimental results (a) and the theo- 
retical calculations for the optimal PIA (b) are depicted 
next to each other. In the theoretical calculations, the 
experimental value is used for the amplitude of the sig- 
nal input. The first and second moments are expressed 
by ellipses, which correspond to the cross sections of the 
Wigner functions. Note that the theoretical ellipses in (b) 
are strictly circles. In each phase space diagram, there 



are three ellipses. The ellipse in green is the signal input. 
Its radius is almost 0.5 which corresponds to the standard 
deviation of vacuum fluctuation. The ellipse in red is the 
signal output. Its center is about y/2 times farther away 
from the origin and its radius is about v3 times larger 
than those of the signal input. The ellipse in blue is the 
idler output. Its radius is almost the same as that of the 
signal output, whereas its center is flipped around the x- 
axes from that of the signal input, which again represents 
the Hermitian conjugate term in Eq. (|2bj) . 

In principle, we can fully characterize our PIA with 
only the above way of analysis that treats quadrature 
values directly. However, such treatment requires a large 
amount of data for good accuracy. Thus, in the follow- 
ing, we resort to power measurements using a spectrum 
analyzer. Not only output quadratures (shown in Figs. [5] 
and [7]) but also their correlations (shown in Figs. |S1 |H 
and II Op are taken for several input states. In each figure, 
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FIG. 6: Output powers for vacuum inputs. Vertical axes are 
powers in dB scale normalized by shot noises. Blue: Output 
quadratures. Cyan: Shot noises. Red: Theory for optimal PIA 
outputs. Magenta: Theory for our PIA outputs with —5 dB 
squeezed ancillas. Green: Theory for our PIA outputs with 
vacuum ancillas. 



results of each quadrature are contained in one of boxes. 
Vertical axes are powers in dB scale which are normalized 
by corresponding shot noises. 

Fig. [6] shows the experimental results for vacuum in- 
puts (fluctuating traces), together with their theoretical 
expectations (straight lines) . There are four boxes corre- 
sponding to the four output quadratures, namely, x° ut , 
p° nt , Xrj ut , an d P2 Ut - The traces in blue are the powers 
of the output quadratures. The traces in cyan around 
dB are the powers of the shot noises, which are used 
for normalization. Since the inputs are in vacuum states, 
the powers of the shot noises correspond to those of the 
input quadratures x 1 / 1 , p™, x™, and p' 2 n . We put three 
kinds of theoretical lines corresponding to three different 
conditions. For the optimal PIA of a coherent state with 
G = 2.0, the output quadrature variances become three 
times larger than the initial shot-noise-limited variance, 
where two from amplification and one from contamina- 
tion by the other mode. The corresponding 4.8 dB is 
marked by the lines in red. Our PIA with finite ancilla 
squeezing is suffered from further excess noise. Assum- 
ing —5 dB of squeezing for the ancilla states, we calculate 
theoretical values which are marked by the lines in ma- 
genta. We also show them with vacuum ancillas by the 
lines in green. The lines in red, magenta, and green are 
very close to each other, thus other experimental errors 
are dominant rather than the ancilla squeezing levels in 
these results. Therefore, it is hard to discuss nonclassi- 
cality of our PIA only from these results. As is shown 
later, the effects of ancilla squeezing more clearly appear 
in the output correlations. 

Next we use several coherent states as inputs. The re- 
sults are shown in Fig. [7] The four input quadratures 
x™, p™, x™, and p™ are displaced from zero mean val- 
ues by turns, leaving the other three quadratures at the 
vacuum level. There are four subfigures labeled from (a) 
to (d) corresponding to such four excitations. For each 
subfigure, there are five boxes. The trace in magenta in 
the leftmost box shows the measured power of the ex- 
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FIG. 7: Output powers for inputs in several coherent states. 
One of four input quadratures x 1 ™, pi", i' 2 n , is excited, at 
the same time the other three quadratures are left at the vac- 
uum level. Vertical axes are powers in dB scale normalized by 
shot noises. Magenta: Excited input quadratures. Red: Out- 
put quadratures with excitation in inputs. Blue: Output 
quadratures without excitation in inputs. Cyan: Shot noises. 
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cited input quadrature. The other four boxes correspond 
to the four output quadratures, namely, i° ut , p° ut , x^ , 
and p2 Ut . The traces in red show the output quadrature 
powers with the input excitation. They are compared to 
those without the excitation shown by the blue traces, 
which arc replottings of the blue traces in Fig. [6] The 
obtained results show the following features. When a 
quadrature x or p of an input mode is excited, the same 
quadratures of both of the output modes are excited, 
whereas the conjugate quadratures do not change from 
the nonexcited levels. The two increased output powers 
differ by about 3.0 dB, where the larger one corresponds 
to the amplified signal and the smaller one corresponds 
to the phase-conjugated idler output. These features are 
exactly what are expected from Eq. ([2]) for G = 2 and 
9 = 0. Note that the coefficients y/2 correspond to the 
3.0 dB. 

The results in Figs. [5] and [7] are only for the five specific 
input states. However, the results for other input states 
can be predicted on the assumption of linearity. More 
precisely speaking, the absolute values of the coefficients 
of i™, p 1 / 1 , x™, and p 2 n in Eq. @ are determined from 
these results. On the other hand, the signs of the coeffi- 
cients are not determined from them. However, they are 
checked from the phase space diagrams shown in Figs. [3] 
and [5] In the above sense, the results shown so far give 
full information of the input-output relation when output 
modes are separately concerned. 

In order to fully characterize our amplifier, the individ- 
ual behaviors of the output modes are not enough. In the 
following, we are concerned with the output correlations. 

Since unitary PIA is equivalent to two-mode squeezing, 
the two output modes should be entangled, and have an 
EPR type of correlation. The results of EPR correlation 
are shown in Fig. [5] Here the two input modes are both 
in vacuum states. There are two boxes corresponding to 
x and p correlations. The lower traces in blue show the 



two-mode squeezing of x° v 



and Px+pl , on the 



other hand, the upper traces in blue show the two-mode 



antisqueezing of x° v 



and Pi — p 2 , respectively. 



They are compared with the summed shot noises of the 
two homodyne detections shown by the traces in cyan. 
Several theoretical lines are plotted together. The lower 
and upper lines in red are the theoretical values of two- 
mode squeezing and antisqueezing for the optimal PIA, 
respectively. Our results of two-mode squeezing are de- 
graded from the ideal case due to finite squeezing of ancil- 
las. Assuming —5 dB squeezing for ancillas, theoretical 
expectation is marked by the lines in magenta. That 
for vacuum ancillas is also marked by the lines in green, 
which is exactly equal to the shot noise level. In contrast, 
the theoretical two-mode antisqueezing is always ideal for 
arbitrary ancillas. The experimental results well agree 
with the theory assuming —5 dB squeezing of ancillas. 
Since the lower traces in blue are both below the traces 
in cyan, existence of entanglement is verified between the 
two output modes via the Duan-Simon criterion [3D . [HJ ■ 
From the nonclassical correlation between the two out- 
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FIG. 8: Two-mode squeezing and antisqueezing for vacuum 
inputs. Vertical axes are powers in dB scale normalized by 
summed shot noises of two homodyne detections. Lower 
Blue: Two-mode squeezing in £? ut - *2 Ut and p° ut + pT* '■ 
Upper Blue: Two-mode antisqueezing in x° ut + X2 Ut and 



f>2 U ■ Cyan: Summed shot noises. Lower Red: The- 



ory of two-mode squeezing for optimal PIA outputs. Upper 
Red: Theory of two-mode antisqueezing for optimal PIA out- 
puts. Magenta: Theory of two-mode squeezing for our PIA 
outputs with —5 dB squeezed ancillas. Green: Theory of two- 
mode squeezing for our PIA outputs with vacuum ancillas. 



put modes, we investigate the reversibility of our PIA. 
For this purpose, we virtually realize the inverse trans- 
formation electrically and reconstruct the initial quadra- 
tures. Neglecting the excess noise from finite ancillas, 
PIA that we demonstrate has the input-output rela- 



tion as a° nt = V2a^ + (af) 



= V2ti£ + (ai n )t, 



which is obtained by substituting G — 2 and = 
into Eq. @. The inverse transformation becomes as 



= V2af - (a 



in ^ t 



V2af - (af)t, 



or equiva- 



lently, 



pT 1 =V2; 



x 2 i 
pt, 



° ut =V2x 2 n 
V2p 2 D 



-out 
V2 



Pi 



(9a) 
(9b) 



Therefore, by adding or subtracting the two homo- 
dyne outcomes with 3.0 dB difference of gains, the ini- 
tial quadratures are reconstructed. The reconstructed 
quadratures are denoted by the superscripts "rec" in the 
following. Note that the initial quantum state is not re- 
covered in the experiment. In addition, note also that 
only one of two quadratures x or p can be reconstructed in 
each moment, and never both simultaneously. The recov- 
ery of the initial state is possible only when either a quan- 
tum channel is between the signal and idler outputs [25[ 
or linear pre-processing is applied [l(|. However, the 
demonstrated reconstruction of the initial quadratures 
can show that necessary correlations for the recovery of 
the initial state are present. 

In Fig. 121 the results of such reconstruction of initial 
vacuum fluctuations are shown. There are four boxes cor- 
responding to the four reconstructed quadratures ii ec , 
Pi° c , i 2 ec , and p r 2 ec . The powers of the reconstructed 
quadratures are shown by the traces in blue. The traces 
in cyan are the powers of the summed shot noises of the 
two homodyne detections, which are taken with the same 
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FIG. 9: Powers of initial vacuum fluctuations reconstructed 
from output correlations. Vertical axes are powers in dB 
scale normalized by shot noises. Blue: Reconstructed quadra- 
tures. Cyan: Summed shot noises of two homodyne detec- 
tions. Red: Theory for optimal PIA. Magenta: Theory for 
our PIA with —5 dB squeezed ancillas. Green: Theory for 
our PIA with vacuum ancillas. 



electric gains as the traces in blue. Note that the lower 
levels of the blue traces than the cyan traces are due to 
the nonclassical correlation. From Eq. the summed 
shot noises should have three times larger variances than 
those corresponding to the initial vacuum fluctuations. 
Thus, we infer the original vacuum level to be 4.8 dB 
below the measured sum of the shot noises. All results 
shown here are normalized by the inferred vacuum level. 
For the optimal PIA, vacuum fluctuations are perfect re- 
constructed, thus the theoretical expectation coincides 
with the vacuum level of dB, which is marked by the 
lines in red. The increases of the blue traces from dB 
show the imperfection of our PIA. The theoretical values 
for —5 dB squeezing of and no squeezing of ancillas are 
marked by the lines in magenta and green, respectively. 
The experimental results of the traces in blue are in good 
agreement with the lines in magenta. 

Next, we pay attention to the reconstruction of mean 
amplitude, using coherent states as inputs. The results 
are shown in Fig. llQI The four input quadratures, namely 
i™, p 1 ", x™ , and p™, are excited one by one, as shown 
in the leftmost boxes of four subfigures. For each excita- 
tion, the powers of the four reconstructed quadratures are 
measured, the results of which are put in the other four 
boxes on the right side. For the leftmost box, the trace 
in magenta is the power of the excited input quadrature, 
and the trace in cyan is the power of the shot noise of the 
corresponding homodyne detection. The increase of the 
magenta trace from the cyan trace indicates the excita- 
tion. For the other four boxes, the traces in red and blue 
show the powers of the reconstructed quadratures with 
and without the input excitation respectively, and the 
traces in cyan are the summed shot noise powers of the 
two homodyne detections. Similarly to Fig. [SJ one third 
of the summed shot noise power is used for normaliza- 
tion. The reconstructed quadratures are excited almost 
to the same levels as the input quadratures at the ex- 
cited quadratures, whereas they remain unchanged from 
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FIG. 10: Powers of initial quadratures reconstructed from 
output correlations. One of four input quadratures k™, pi", 
$2 ) P2 1 is excited, at the same time the other three quadra- 
tures are left at the vacuum level. Vertical axes are powers in 
dB scale normalized by shot noises. Magenta: Excited input 
quadratures. Red: Reconstructed quadratures with excitation 
in inputs. Blue: Reconstructed quadratures without excita- 
tion in inputs. Cyan: Shot noises of a single homodyne detec- 
tion for the input quadratures and two homodyne detections 
for the reconstructed quadratures. 
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the nonexcited levels at the nonexcited quadratures. 

All the results shown above prove the success of our 
demonstration of PIA. They agree well with the theoret- 
ical calculations assuming —5 dB of squeezing for ancillas. 



VI. EXPERIMENTAL RESULTS FOR CLONER 

We show next the results of the cloning experiment in 
a manner similar to the PIA experiment, i.e., quadrature 
data and the phase space diagrams reconstructed from 
them are used for intuitive understanding and also for 
check of the mirror image that is found in the anticlone, 
and the full verification is given by the power analysis 
for various input states. In the following, since the rep- 
resentation and interpretation of results are almost the 
same as those for the PIA experiment, we give only short 
description of them. The signal input of PIA is denoted 
by the term "original" , and the resulting two clones and 
one anticlone are denoted by "clone-1", "clone-2", and 
"anticlone" , respectively. There are also abbreviation of 
these terms as "org", "cln-1", "cln-2", and "a-cln" in the 
figures and mathematical expressions. Occasionally, the 
term "input" is used indicating the original, and "out- 
put" indicating the two clones and one anticlone. 

Fig. [TT] shows the quadrature data. From an origi- 
nal (a) in a coherent state, two clones (b) and (c) and an 
anticlone (d) are produced. We see that the original and 
the two clones have almost the same sinusoidal curves of 
the mean amplitudes, though the fluctuations are uni- 
formly increased in the clones. On the other hand, the 
anticlone has the same sinusoid when the phase is flipped. 

Fig. [12] is the phase space diagram computed from the 
quadrature data shown in Fig. [TT] The first and second 
moments of each distribution are represented by an el- 
lipse. Next to the experimental diagram (a), theoretical 
calculations for the optimal cloning is depicted (b). The 
two ellipses of the clones (red and magenta) in the ex- 
perimental diagram are almost overlapped. The centers 
of the two clones are almost the same as that of the orig- 
inal (green), whereas the radii of the clones are larger 
than that of the original. On the other hand, the anti- 
clone (blue) has a different center, where the sign of p 
quadrature is opposite to that of the original. 

We move on to the power analysis. First we show the 
powers of the output quadratures (shown in Figs. [TBI and 
[TT]) . and then their correlations (shown in Figs. [T5l [TBI 
and [17]) . Fig. [T3l shows the cloning of a vacuum state, 
and Fig. [TT] shows the cloning of several coherent states. 

In Fig. 1131 there are six boxes corresponding to six 
output quadratures, namely, x c i n -i, p c in-i, icin-2, Pdn-2, 
^a-cin , and Pa-cin ■ For each box, there are two experimen- 
tal traces. The traces in blue are the powers of the output 
quadratures. The traces in cyan are the shot noise pow- 
ers used for normalization, which are equal to the powers 
of the input quadratures. There are also three kinds of 
theoretical lines. The lines in red are for the optimal 
cloning, and the lines in magenta and green are for our 
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FIG. 11: Quadrature data. The phases of homodyne mea- 
surements are scanned from to 2tt, which correspond to 
horizontal axes. Vertical axes are quadrature values. 
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(a)Experiment. (b)Theory. 

FIG. 12: Phase space distributions computed from the 
quadrature data in Fig. [TTJ The first and second mo- 
ments of Gaussian Wigner functions are expressed by ellipses. 
Green: Original. Red: Clone-1. Magenta: Clone-2. Blue: An- 
ticlone. 



cloning using —5 dB squeezed and vacuum ancillas for 
PIA, respectively. Note that the lines in red at 3.0 dB 
for clones correspond to the cloning limit. From these re- 
sults, cloning fidelities are estimated at F = 0.63 ± 0.01 
for a vacuum original, which is higher enough than the 
classical limit of F — 1/2 and very close to the cloning 
limit of F = 2/3. 

In Fig. [TTJ there are two subfigures (a) and (b) cor- 
responding to excitations in x org and j5 rg, respectively. 
Each subfigure is composed of seven boxes, where the 
leftmost one shows the input excitation in i org or p OTg , 
and the other six boxes show the output quadratures of 
icin-i, Pcin-i, icin-2, Pcin-2, ^a-cin, and p a -cin- The trace 
in magenta shows the power of the excited input quadra- 
ture. The traces in red and blue show the output powers 
with and without the input excitation, respectively. The 
traces in cyan are the shot noise powers used for normal- 
ization. When we excite one quadrature of the original, 
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FIG. 13: Output powers for vacuum originals. Vertical axes 
are powers in dB scale normalized by shot noises. Blue: Out- 
put quadratures. Cyan: Shot noises. Red: Theory for optimal 
doner. Magenta: Theory for our doner with —5 dB squeezed 
ancillas. Green: Theory for our doner with vacuum ancillas. 



(a)EPR correlation between clone-1 and anticlonc. 
x P 



1 



CD 



16 



12 



3 
o 

Q_ 



x org x cln-1 Pcln-1 x cln-2 Pcln-2 x a-cln Pa-cln 



16 



12 



























































































Porg 


(a)(£ org ) ^ 0. 
x cln-1 Pcln-1 x cln-2 Pcln-2 


x a-cln 


Pa-cln 



S -1 
s 

o 

o_ -2 
-3 
-4 




(b)EPR correlation between clone-2 and anticlone. 

FIG. 15: EPR correlation between each clone and the 
anticlone. Vertical axes are powers in dB scale normal- 
ized by summed shot noises of two homodyne detections. 
Blue: Xdn — ia-cin and p c in + p a _ c i n . Here the subscripts "cln" 
denote clone-1 for (a) and clone-2 for (b). Cyan: Summed shot 
noises of two homodyne detections. Red: Theory for optimal 
doner. Magenta: Theory for our doner with —5 dB squeezed 
ancillas. Green: Theory for our doner with vacuum ancillas. 



and that between the clone-2 and the anticlone is shown 
in Fig. 



15(b) 



By electrically adding or subtracting 
the two homodyne signals with the same electric gains, 
four observables are measured, namely, i c in-i — ia-cin, 

Pcln-1 + Pa-cln, ^cln-2 ~ ia-cln, and p cln _ 2 + pa-cln, which 

are separately contained in boxes and shown as the blue 
traces. These traces are all below the summed shot noises 
of the traces in cyan. From these results, we verify bipar- 
tite entanglement between each clone and the anticlone 
from Duan-Simon criterion [3ll l32|. and eventually, tri- 
partite entanglement of Class I where none of three par- 
tial systems is separable from the others @ . Theoretical 
lines in red, magenta, and green are plotted together, 
corresponding to infinite squeezing, finite squeezing of 
—5 dB, and no squeezing of ancillas for PIA, respectively. 
The experimental results of the traces in blue well agree 
with the lines in magenta. 

Using the nonclassical correlations, we reconstruct the 
original quadratures. The reconstructed quadratures are 
denoted by a; rcc and p lcc . The results for a vacuum state 
are shown in Fig. [T6l and those for coherent states are 
shown in Fig. 1171 For the reconstruction, three homo- 
dyne signals are added with the same electric gains and 
tween the clone-1 and the anticlone is shown in Fig. 15(a) appropriate signs. For the similar reason to the PIA ex- 



(b){p org > #0. 

FIG. 14: Output powers for originals in several coherent 
states. One of two original quadratures ^org, Porg is excited, at 
the same time the other quadrature is left at the vacuum level. 
Vertical axes are powers in dB scale normalized by shot noises. 
Magenta: Excited original quadratures. Red: Output quadra- 
tures with excitation in originals. Blue: Output quadratures 
without excitation in originals. Cyan: Shot noises. 



the same quadratures in the three output modes are ex- 
cited to almost the same level, whereas the conjugate 
quadratures do not change. 

Our remaining concerns are the output correlations 
and the reversibility. First, in Fig. [151 we show EPR 
correlation between each clone and the anticlone as a 
sufficient condition for entanglement. The correlation be- 
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FIG. 16: Powers of initial vacuum fluctuations reconstructed 
from output correlations. Vertical axes are powers in dB scale 
normalized by shot noises. Blue: Reconstructed quadratures. 
Cyan: Summed shot noises of three homodyne detections. 
Red: Theory for optimal doner. Magenta: Theory for our 
doner with —5 dB squeezed ancillas. Green: Theory for our 
doner with vacuum ancillas. 
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noise power is used for normalization. 

In Fig. [T6l the powers of the reconstructed vacuum 
fluctuations are plotted as the traces in blue, and com- 
pared to that of the summed shot noise plotted as the 
traces in cyan. The blue traces are below the cyan traces 
for both x Tec and p rec due to nonclassical correlations. 
Theoretical expectations are also shown as the lines in 
red, magenta, and green corresponding to the three dif- 
ferent conditions of infinite squeezing, finite squeezing of 
—5 dB, and no squeezing of ancillas, respectively. The 
perfect reconstruction corresponding to dB is marked 
by the red lines, which is not achieved in the experiment 
due to the finite squeezing of ancillas. The results are 
degraded almost to the level of the magenta lines, as ex- 
pected from the theory. However, they are still slightly 
below 3.0 dB which corresponds to the cloning limit. 
From these results, the fidelity of reconstruction is calcu- 
lated for a vacuum state. A perfect unitary cloning allows 
the reconstruction fidelity of F — 1. The experimental 
value is calculated as F — 0.74 ±0.01, which is higher 
than the cloning limit of F = 2/3. The cloning limit can 
be considered as the classical limit for the reproduction 
of the original state, because one can never obtain a bet- 
ter approximation of the original state than the clones 
if the nonclassical correlations between the clones and 
anticlones can not be utilized. 

In Fig. [TTJ the input quadratures x OIg and p OIg are ex- 
cited one by one, and for each case they are reconstructed 
from the output correlations. The trace in magenta in the 
leftmost box is the power of the excited original quadra- 
ture, whereas the trace in cyan in the same box is the 
power of the corresponding shot noise. The traces in red 
and blue in the other two boxes on the right side are the 
powers of the reconstructed quadratures with and with- 
out the excitation, respectively, and the traces in cyan at 
4.8 dB are the powers of the summed shot noises of the 
three homodyne detections. At the excited quadrature, 
almost the same level of excitation is reconstructed. In 
contrast, at the conjugate quadrature, no effect of the 
excitation is observed. 



(b)(p org > #0. 



VII. SUMMARY 



FIG. 17: Powers of initial quadratures reconstructed from 
output correlations. One of two original quadratures x r S , Porg 
is excited, at the same time the other quadrature is left at the 
vacuum level. Vertical axes are powers in dB scale normal- 
ized by shot noises. Magenta: Excited original quadratures. 
Red: Reconstructed quadratures with excitation in originals. 
Blue: Reconstructed quadratures without excitation in orig- 
inals. Cyan: Shot noises of a single homodyne detection for 
the input quadratures and three homodyne detections for the 
reconstructed quadratures. 



periment, the summed shot noise has three times larger 
variance than that corresponding to the vacuum fluctua- 
tion of the original. Thus, one third of the summed shot 



We succeeded in phase-insensitive optical amplification 
in a reversible manner. Our amplifier preserves the idler 
output. The entanglement between the signal and idler 
is responsible for the reversibility. The scheme is ba- 
sically based on linear optics, homodyne measurements 
and feedforward. Offline-prepared squeezed states which 
are used as ancillas provide nonclassical properties for 
our PIA. We demonstrated for the amplification gain of 
G = 2.0. By splitting the amplified output in half, we 
also demonstrated 1 — > 2 approximate cloning of coherent 
states, where the remaining idler output was interpreted 
as the anticlone. 

For both experiments, the full demonstration was given 
in the following sequence. First, we characterized the in- 
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dividual output modes. By treating the quadrature data 
directly, they were visualized as phase space diagrams to 
help intuitive understanding. Especially, the mirror im- 
age in the idler output or the anticlone is shown. Then 
the input-output relation was examined more strictly by 
using several different coherent states as inputs. Finally, 
the output correlations were examined. They are impor- 
tant because the nonclassical properties are only acces- 
sible via them. Not only the ordinary EPR correlations 
were shown, but also the possibility of the reverse opera- 
tion was directly presented by appropriate measurements 
of the correlations. 

Our results are a good demonstration that shows the 
properties of an amplification process, which have been 
theoretically known for decades but not fully demon- 
strated experimentally. Especially, it is reversible when 
the idler is present. Such reversible amplification is signif- 
icant from practical respects, as is shown in Refs. [Tol l25j|. 
We did not demonstrate the reverse operation but showed 
its possibility from the correlations. The recovery of the 
signal state only requires the Bell measurement and feed- 
forward [lol [23, which would be much less lossy and 
noisy than implementing the inverse transformation. Full 
and partial recovery of the distributed information via 
such a feedforward scheme is left for future experiments. 
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Appendix A: General number of clones 

In this appendix, the discussion in Sec. Illll is extended 
to K — ► L cloning. 

The procedure of K — > L symmetric cloning can be de- 
composed into three steps as follows J2J [22| . Firstly, all 
the information of (aidjPd) & r e put together into a single 
mode by a beamsplitter network. A state with larger 
amplitude D(^/Kxd, VKpd)\^P) is created from the K 
identical originals D(xd,pa)\ip) in this step. Secondly, 
the combined signal is amplified with the gain G = L/K. 
Finally, the amplified signal is combined with L — 1 ancil- 
las by another beamsplitter network, creating L clones. 
For asymmetric cloning, the procedure is essentially the 
same, but the amplification in the second step is applied 



to a part of the combined signal, and the gain is changed 
correspondingly as G = 1 + J2k=i Uk B3' ( - )n tne other 
hand, L — K anticlones are obtained from the idler out- 
put by combining it with L — K — 1 ancillas by another 
beamsplitter network. Therefore, as a whole, L clones 
and L — K anticlones are obtained, whose complex mean 
amplitudes are a and a* respectively, from K originals 
with the complex mean amplitude of a. It is evident 
that there is no entanglement among clones or anticlones, 
however, there is entanglement between a clone and an 
anticlone. 

The asymmetric clones have the following form: 

L-l 

ficln-fc = "T^Oorg + V^^idl + J! K fc^anc-£, (Al) 

k=l 

where a c i n -fe is the annihilation operator of the fc-th clone, 
a' is that of the combined original after the first step, 
fiidi is that of the idler input, and a anc _i is that of the £-th 
ancilla input for the latter beamsplitter network. The 
coefficients y^nfc and Kki are not independent in order 
to preserve the commutation relations. By setting all 
the ancillas in vacuum states, the added noises become 
rotationally symmetric and Gaussian, and their variances 
correspond to the parameters n k [33J. They satisfy the 
relation below 16]: 

(EV^) 2 = (i- K) (jr n k + l). (A2) 
fc=i fe=i 

The optimality of Eq. (|A2[) is proven with respect to the 
cost function constructed from the variances |16j : 

L 

C(rii , . . . , n L ) = c k nk- (A3) 

k=l 

In particular, for the symmetric cloning, i.e., m = • • • = 
til = n, Eq. (|A2[) saturates the following inequality which 
is obtained in Ref. |34] from the consistency with the 
uncertainty relation: 

( A4 ) 

The limit fidelity for symmetric Gaussian cloning is 
calculated as F = KL/(KL — K + L). On the other hand, 
taking the limit that L goes to infinity, the classical limit 
of cloning (i.e., the limit of state estimation) is obtained 
as F = K/(K + l). 

As with 1 — > 2 cloning, optimal 1 — > L cloning can 
be fully reversed by L — 1 Bell measurements performed 
on each set of a clone and an anticlone and subsequent 
feedforward to the remaining single clone [25[ . Even with 
the smaller number of Bell measurements, the original is 
partially recovered accordingly. 
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